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In this paper we study the relation between long cycles and Bose-Condensation in the 
Infinite range Bose-Hubbard Model with a hard core interaction. We calculate the 
density of particles on long cycles in the thermodynamic limit and find that the exis- 
tence of a non-zero long cycle density coincides with the occurrence of Bose-Einstein 
j — , condensation but this density is not equal to that of the Bose condensate. 
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1 Introduction 



In 1953, Feynman analysed the partition function of an interacting Bose gas in terms of 
the statistical distribution of permutation cycles of particles and emphasized the roles of 
long cycles at the transition point [7j. Then Penrose and Onsager, pursuing Feynman's 
arguments, observed that there should be Bose condensation when the fraction of the to- 
tal particle number belonging to long cycles is strictly positive |12j . These ideas are now 
commonly accepted and also discussed in various contexts in systems showing analogous 
phase transitions such as percolation, gelation and polymerization (see e.g. [I], [T^j, [E]), 
though it has been recently argued by Ueltschi [19] that in fact the hypothesis is not always 
valid. To our knowledge, there had not appeared a precise mathematical and quantitative 
formulation of the relation between Bose condensate and long cycles until the work of Siito 
[To] and its validity has been checked only in a few models: the free and mean field Bose gas 
in [16], (see also Ueltschi [18]) and the perturbed mean field model of a Bose gas studied in 
[5]. In these models it is shown that the density of particles in long cycles is equal to the 
Bose condensate density. The purpose of this paper is test the validity of the hypothesis in 
yet another model of a Bose gas, the Infinite range Bose-Hubbard Model with a hard core. 
Here we calculate the density of particles on long cycles in the thermodynamic limit and 
find that though the existence of a non-zero long cycle density coincides with the occurence 
of Bose-Einstein condensation, this density is not equal to the Bose condensate density. 

The main simplifying feature in this model is the following. In general the density of particles 
on cycles of length q for n particles can be expressed (apart from normalization) as the trace 
(see for example Proposition 13.11) of the exponential of the Hamiltonian for n — q bosons 
and q distinguishable particles (no statistics). In terms of the random walk representation 
(cf p2]), the particles in this model are allowed to hop from one site to another with equal 
probability. We can prove (Proposition 13.21) that in the thermodynamic limit we can neglect 
the hopping of the q particles so that bosons have to avoid each other and the fixed positions 
of the distinguishable particles. This is equivalent to a reduction of the lattice by q sites. 
Moreover the q particles are on a cycle of length q. For q > 1, this means for example, that 
the position of the second particle at the beginning of its path is same as the position of 
the first particle at the end of its path. But since they do not hop this is impossible by 
the hard core condition and therefore among the short cycles only the cycle of unit length 
contributes. Since the sum of all the cycle densities gives the particle density, this means 
that in the thermodynamic limit the sum of the long cycle densities is the particle density 
less the one-cycle contribution. The one-cycle density, apart from some scaling and the 
normalization, is then the partition function for the boson system with one site removed 
from the lattice, which can be calculated. 

The model without a hard-core will be treated in another paper. There we can again neglect 
the hopping of the q distinguishable particles. However in that case cycles of all lengths 
contribute to the long-cycle density. It is relatively easy to see that when there is no con- 
densation the long-cycle density vanishes but we do not yet know what happens when there 
is Bose-Einstein condensation. 

In Section 2 we first describe the model and recall its thermodynamic properties as stated by 
Penrose [13] (see also Toth [17] and Kirson[9]). We then apply the general framework for cycle 
statistics described in [3], following [TTJ. Using standard properties of the decomposition of 
permutations into cycles, the canonical sum is converted into a sum on cycle lengths. This 



makes it possible to decompose the total density p = p s hort + Piong into the number density 
of particles belonging to cycles of finite length (p s hort) an d to infinitely long cycles (piong) in 
the thermodynamic limit. It is conjectured that when there is Bose condensation, pi ong is 
different from zero and identical to the condensate density. The main purpose of the paper is 
to check the validity of this conjecture in our model. At the end of Section 2 we state in the 
main theorem describing the relation between Bose-Einstein condensation and the density 
of long cycles for our model. 

In Section 3 we prove the main theorem and in Section 4 we discuss briefly Off- diagonal 
Long-Range Order. 



2 The Model and Results 

The Bose-Hubbard Hamiltonian is given by 

# BH = j K-<)K-«J + A E n >-- 1 ) ( 21 ) 

x,y£Ay :\x— J/|=1 x&Ay 

where A v is a lattice of V sites, a* and a x are the Bose creation and annihilation operators 
satisfying the usual commutation relations [a*, a y ] = 5 x>y and n x = a* x a x . The first term with 
J > is the kinetic energy operator and the second term with A > describes a repulsive 
interaction, as it discourages the presence of more than one particle at each site. This model 
was originally introduced by Fisher et al. [8]. 

The infinite-range hopping model is given by the Hamiltonian 

H ™ = ^7 E « " a W a * ~ a v) + A E n *( n * - !)■ ( 2 - 2 ) 

x,y&kv xeAy 

This is in fact a mean-field version of (12. ip but in terms of the kinetic energy rather than the 
interaction. In particular as in all mean-field models, the lattice structure is irrelevant and 
there is no dependence on dimensionality, so we can take A v = {1, 2, 3, ... , V}. The non- 
zero temperature properties of this model have been studied by Bru and Dorlas [3] and by 
Adams and Dorlas [TJ. We shall study a special case of (12. 2p . introduced by Toth [T7j where 
A = +oo, that is complete single-site exclusion (hard-core). The properties of this model 
in the canonical ensemble were first obtained by Toth using probabilistic methods. Later 
Penrose [13] and Kirson [9] obtained equivalent results. In the grand-canonical ensemble the 
model is equivalent to the strong-coupling BCS model (see for example Angelescu [2]). Here 
we recall the thermodynamic properties of the model in the canonical ensemble as given by 
Penrose. 

For p e (0,1), let 



In '- if p^ 1/2, 

if p = 1/2. 

For each (3 > 2 the equation (3 = g(p) has a unique solution in (0, 1/2] denoted by pp (see 
FigHJ). We define p p := 1/2 for (3 < 2. 



g(p)k 




Figure 1: Definition of p fj 
Theorem 2.1 (Penrose [10], Theorem 1) 

The free energy per site at inverse temperature (3 as a function of the particle density p £ 
[0, 1], fp(p), is given by 

(p+-(plnp+(l-p)ln(l-p)) if pe [0, P p]U [1-^,1], 

fM = { 2 p i 

yp + p p {l- p p ) + -{pp\np p + {l- pp)\n{\- pp)) if p e [p(3,i - pp\. 

The density of particles in the ground state in the thermodynamic limit is given by 



n/v=p z,yeA v 



(2.3) 



where ( • ) denotes the canonical expectation for n particles. Penrose showed that for certain 
values of p and (3, Bose- Einstein condensation occurs, that is, p c p > 0. The Bose-condensate 
density is given in the following theorem. 

Theorem 2.2 (Penrose [10], Theorem 2) 

The Bose-condensate density, p c p at inverse temperature [3 as a function of the particle density 
p £ [0, 1], is given by 



if p£ [0, P/3 ]U[l-p^,l], 

{p-pp){l -p-pp) if pe\pp,l-p p \. 



We note that both fp(p) — p and the condensate density pp are symmetric about p = 1/2. 
This can easily seen by interchanging particles and holes. The Boson states being symmetric 
can be labelled unambiguously by the sites they occupy but equivalently they can be labelled 
by the sites they do not occupy (holes). 



Before proceeding to the study of cycle statistics we need to define the n-particle Hamiltonian 
more carefully. The single particle Hilbert space is Ti, v := C v and on it we define the operator 

Hy=I- Py 

where P v is the orthogonal projection onto the unit vector 

gv = -=(l,l,...,l)eW v . 
W 



In terms of the usual basis vectors of H v , {ej | i — 1 . . . V}, P v is given by 

p _lf 

i=i 

Thus H v is the orthogonal projection onto the subspace orthogonal to g v . For an operator 
A on H v , we define A (n) on = H V ®H V ®---® H Vj by 

n times 

A {n) = A® I ® . . .® I + 1 ®> A® . . .® I + . . . + 1 ® I ® . . .® A. 

(n) 

With this notation we can define the non-interacting n-particle Hamiltonian H v acting on 
the unsymmetrised Hilbert space T~Cy ^ as: 

H (n) = l{n) _ p (n) 

= n-Py®I®---®I-I®Py®I®---<g>I I ® I ® ■ ■ ■ ® Py. 

For bosons we have to consider the symmetric subspace of TCy . The symmetrisation pro- 
jection cr™ on TCy is defined by 



where U n : Hy l— > "Hv is the unitary representation of the permutation group S n on Hy 
defined by 

f^rOl ® 02 ® • • • ® 0») = 0tt(1) ® 0tt(2) ® • • • ® 07r(n), 0j £ W v , j = 1, • • • , U\ 7T G S^. 

Then the symmetric n-particle subspace is := a™Hy . 
When Hy^ is restricted to Tly + , we obtain 

( a l~ a *y)( a x- a y)- 

We introduce the hard-core interaction by applying a projection to Tly to forbid more than 
one particle from occupying each site. Let {ei}J =l be the usual orthonormal basis for H v . 
We then define the hard core projection on TCy by 

^ hr/ J if ej, = e,:. for some k ^ k' 

7>£ c (e il ®e ia ®..-®e in ) = S „ ^ + u ■ ( 2 - 5 ) 

• • ■ ® e iri otherwise. 



We shall call Hl c 



V^H^ the unsymmetrised hard-core n-particle space and 

P^ c 7i^+ the symmetric hard-core n-particle space. Note that as [U n , V^ ] = 

(J,ri r 



for all 7r G S n , commutes with the symmetrisation and so 7i^ v + 
The hard-core n-particle Hamiltonian is then 

rrhc T>hc trW-phc 

n n,v ■— 1 n 11 V I n 



(2.6) 



acting on the hard-core n-particle space Ti.^ c v . Therefore the Hamiltonian for the infinite- 
range Bose-Hubbard model with hard-core is (12. 6p acting on the symmetric hard-core 
n-particle space Ti^ v ,. 

We shall now analyse the cycle statistics of this model. 

Using (12. 4p . the canonical partition function for the hard-core boson model may be written 

as 



Z p (n,V) = trace H h 



n,V,+ 



-/3H hC ,r 

P " n, v 



n n,V 



1 



— > trace -xyhc 

f)) Z— / L n,\ 



TT£S n 



Following [5] , we define a probability measure on the permutation group S n by 



P"(tt1 



Z (n,V)n\ 



(2.7) 



From the random walk formulation (see for example [T7j) one can see that the kernel of 
Q-p n ,v i s positive and therefore the righthand side of (12. 7p is positive. 

Each permutation it G S n can be decomposed uniquely into a number of cyclic permutations 
of lengths qi, q 2 , ■ ■ ■ , q r with r < n and q\ + q 2 + • • • + q r = n. For q G {1,2,..., n}, let N q (ir) 
be the random variable corresponding to the number of cycles of length q in ir. Then the 
expectation of the number of g-cycles in the canonical ensemble is: 

r=l 

and the average density of particles in g-cycles for the system of n bosons is 

qE v (N q ) 



This brings us then to the following definition. 



V 



(2.8) 



Definition 1 The expected density of particles on cycles of infinite length is given by 

n 

= lim lim £ c v (q). (2.9) 

n/v=p q=Q+l 



For the free Bose gas, the mean field and the perturbed mean field Bose gas, it has been 
shown that pjj° ns = p° the condensate density. For our model, the situation is different. 
Below we state the main result of this paper: 



Theorem 2.3 The expected density of particles on cycles of infinite length, p° ng , at inverse 
temperature (3 as a function of the particle density p £ [0, 1], is given by 



We note that (see Fig 12]): 

• p l ° ng = if and only if p c p = 0. 

• p l ° ng is not symmetric with respect to p = 1/2. As mentioned above the symmetry 
of the model about p = 1/2 is due to the particle- hole symmetry. But the equivalent 
labelling of states by sets of occupied or unoccupied sites (particles and holes) cannot 
be used for distinguishable particles. We shall see (Proposition 13.11) that the g-cycle 
occupation density c"(g) involves q distinguishable particles and n — q bosons and 
therefore the particle-hole symmetry is broken. 

• When p c > 0, p l ° ng starts below p c p since its slope at pp is equal to 1 — 2p while p° p has 
slope 1 — (3p and (3 > 2. Conversely, p l ° n& finishes above p c fj since its slope at 1 — p is 
less than that of p\. 





Figure 2 



and p ° ng for (3 > 2 



3 Proof of the Main Result 



In this section we shall prove Theorem 12.31 First we note that if n/V = p, then 




(3.1) 



so that if we define 



we have 



p s 3 hort = lim lim V c£(g) 

Q— >0O n,V-»oo ' * 

n/v=p 9=1 



long _ short 



For p^ hort we can interchange the sum over q and lim n ,v^oo, 

n/v=p 



pf> rt = lim £ lim 

9=1 n/v=p 

making it much easier to calculate. In fact we shall prove that: 



jp if [0,p„]U [1-^,1], 

if P e[ P0 ,l-p,}. [6 - 2) 



The proof is in four steps. The first step is to obtain a convenient expression for c™(q), 
the mean density of particles belonging to a cycle of length q. We shall denote the unitary 
representation of a g-cycle by U q : Hv^ — ► "Hv\ that is 

Uq^n ® ■ • ■ ® <pi n ) = <pi 2 ® • • • 0i 9 ® 0n- 

When there is no ambiguity we shall use the same notation U q for t/ g ® : ft^ -> 7Y^ n) 

where / is the identity operator. Note that [t^,^ ] = and [t/ ? ,cr"] = 0. 



Proposition 3.1 



c n v {q) 



Z p (n, V)V 



trace n 



he 



U q e 



where H^ n v 



(9) 

V 



n 



Note that though we write this proposition for our special case, in fact c™(g) can be expressed 
in this form for any Boson model with a symmetric Hamiltonian. 

By using cycle statistics, we split our symmetric hard-core Hilbert space TC^ V + into a tensor 

product of two spaces, an unsymmetrised g-particle space Tiy and a symmetric n — q particle 
space Hy+ q \ with the hard-core projection applied. Writing 

A® := A iq) <8> I (n ~ q) and A in ~ q) := I (q) ® A {n ~ q) 

for any operator A on 7i v , we can express our Hamiltonian (12.61) on 7Yg C n v as follows: 

fljv = -Pn (n - 4 9) - ^ 
Let = V^Py and define the following reduced Hamiltonian 

ffj^v = Pn (n - PM) (3.3) 



so that 



rrhc _ rrhc p(<?) 

n,V 11 q,n,v r V • 



-j3H hc 



X,n,V 



The next step is to estimate the effect of neglecting the action of the P v q ^ term (equivalent 
to the hopping of the q particles) in the unsymmetrised space. Let 

c v (q) = ^l^Itrace^ [u q e~^] , (3.4) 

and define 

Z (\, n, V) = trace W hc^ 
where 

A&v = ^° (n - AP< n) ) (3.5) 
Then we have the following estimate. 

Proposition 3.2 

ic^)-^)i< (1 -;' ft) ^ ( ^"-^- g) 

1/ Zp{n,V) 

In the third step we obtain the limit of the ratio on the righthand side of the last inequality: 
Proposition 3.3 

Hm Z p (^,n-q,V-q) _ ( ffl e*> if p G [0, p ] U [1 - p fj , 1], 

The final step is a simple proposition where we check the following: 
Proposition 3.4 Cy(q) — if q > 1 and 

hm c«(l) = { P B( , if ^[0, Pi3 ]U[l-p„l], 

Using these four results the main result, Theorem 12.31 follows very easily. From Propositions 
13.21 and 13731 we have 

lim Cy(q)= lim c v (q). 

n,V — ►oo n,V — >oo 

n/v=p n/v=p 

Since by Proposition 13.41 c"(g)=0 if q > 1, it follows that 



pf rt = lim c »(l] 

n, V — »oo 

n/v=p 



p if p G [0,^] U [1-08,1], 

P ,S"-^ if p g [p, fl , 1 - p,] 



which is the required result. 

In the next four subsections we prove the results stated above. 



3.1 Proof of Proposition I3TTI 



We recall the following facts on the permutation group. 



The decomposition into cycles leads to a partition of S n into equivalence classes of per- 
mutations with the same cycle structure C q , where q = [qi,q 2 , . . . , q r ] is an unordered 
r-tuple of natural numbers with qi + q 2 + ■ ■ ■ + q r = n. 

Two permutation tt' and tt" belong to the same class if and only if they are conjugate 
in S n , i.e. if there exists a tt E S n such that 



// -l / 

TT = TT TT TT. 



(3i 



The number of permutations belonging to the class C q is 



n q !(gig 2 - -.q r 



(3.9) 



with n q ! = ni\n 2 \ . . . rij\ . 



and rij is the number of cycles of length j in q. 



We observe that since our Hamiltonian is symmetric {[H^ v , U n ] = 0,7r E S n ) and therefore 
for tt', tt" E C q , one has by 



trace T^hc 



tr£LCG T-/hc 

ri n,V 



(3.10) 



For q E N, let N q (7r) be the number of cycles of length q in tt. 

Let Tj denote the number of cycles of length j. Then Y2jj r j = n anc ^ the corresponding 
number of permutations this cycle structure is n\/ Y\j j Tj rj\ (from ( 13. 9ft ). Denote (rj) the 
class of permutations with such a cycle structure. Then 



K(N q = r) 



Z p (n, V)n\ 



Yl trace *£ v 



(rj) 7!-e(rj) 
r„=r 



Z f (n,V)n\^U,>J T 



E 



ill 



~tr£LCG ay 



(rj) 1Ll J ' 



he 

n,V 



n,V 



where tt is any permutation with cycle distribution (rj). Suppose that r > 1 and consider 
a permutation where the first q indices belong to the same cycle of length q. Let tt' denote 
the permutation of the remaining n — q indices. We have 



U w = U q ® U w > 



and 7r' has cycle structure (rj — 5j q ). Then 

1 1 ^ (n-q)\ 



K(N g = r) 



Z (n, V)n\ 



/ ~v~t '■ i 1 rrtrace 

[[ j>1 3 rj rj\ qr(n — qy. n < v 



(U q ® t^)e-^»> 



r a =r— 1 



E 



Zp(n,V)qr(n-q)\ 

7T tOn-g 



Then the canonical expectation of the number of g-cycles is found to be 

oo 



r=0 



Zp(n,V)q(n-q)\ , 



E 



n'es„ 



7( l T/ , ( 1 r. £ trace wW ® ^)e-^-^ c 

7T kzOn — q 



1 1 

-trace 



1 1 



Z p (n, V) g tr ace v% (H^ } ) 



„- 9) te c (f/ 9 ®/ (n - 9) )e-^V^ c 
L 

(17, ® J^V"^ 



Since 



we have proved Proposition 13. II 



3.2 Proof of Proposition 15721 



qK(N q 
V 



To prove Proposition 13.21 we have to obtain an upper bound for 

■ trace 



trace oyhc 

q,n,v 



U q e~^v 



{ hc 

1,n,V 



In order to do this we first shall introduce some notation and make some remarks before 
proceeding. 

Let A[r + q ^ be the family of sets of n — q distinct points of Ay. For k = {k\, /c 2 , . . . , k n _ q } G 
A^- q) let 



a+ "(e kl ® e fc2 ® • • ■ ® e fcn _ 

/he 

L (n-(?),v,+ ' — r n-q' 1 -V+ 



Then {|k> | k e A?+ q] } is an orthonormal basis for HY_ „, „ , := Vll n H^~ q) 



Similarly let Ay be the set of ordered g-tuples of distinct indices of A v and for i 
(ix,i 2 , ...,i q )e A^ let 



Then {|i) | i 6 Ay } is an orthonormal basis for H^ v := V\ C U^' . 

If k G Ay+ and i G Ay we shall write k ~ i if {k±, /c 2 , • • • , &n- g } H {z'i, Z2, • • • , z 9 } = and 
we shall use the notation 

|i;k) := |i) g> |k). 

Then a basis for Wj^ n v may be formed by taking the tensor product of the bases of Ti.^_ q ^ v + 
and 7~L\ C V where we disallow particles from appearing in both spaces simultaneously. Thus 
the set {|i; k) | k G A<?~ q \ i G A$\ k ~ i} is an orthonormal basis for Ti.„„ ,,. 

We shall need also the following facts. For simplicity we shall write H and P for H q c n v and 
P^ 1 respectively, as defined in equation (13. 3ft . 

Let Pj" ^ be the projection of 7i^_ q -, v + onto a space with none of the n — q particles at 
the points z'i, %% ■ ■ ■ , i q (so there are V — q available sites for n — q particles) and not more 
than one particle at any site. Then 

Remark 3.1 For i ~ k, if s > 

e- pS '\i; k) = |i; e^'tye - ^* (3.11) 
where H' 1 = v[ n ~ q \{n - q) - P^'^P^. 

This can be seen as follows: For i ~ k, 

tf|i;k>=^(n-P^>n C |i;k> 

= qV^\i ] k)+V^\r,((n-q)-P v n ~ q) )k) 
= <z|i; k) + |i; V[ n - q \{n - q) - P v n ~ q) )k) 
= q\i; k) + \i;H i k). 



Remark 3.2 For i ~ k, 



1 n-q 

H l \k) = (n - q)\k) - - K*i' k *> ■ ■ ■ ^ • • • ' ( 3 - 12 ) 



i=l J^iUk\{fc,} 



where the hat symbol implies the term is removed from the sequence, while from A3.5\) . for 
k G HyZ q + we have 



H 



he 

X,n—q,V—q 



A 



n—q V—q 



k) = (n-q)\k)--—J2 E \(ki,k 2 ,...Xj,---,k n - q )). (3.13) 



H 1 = 1 3=1 



Thus H 1 is unitarily equivalent to H^_ q y v n _ q v _ q and 



trace oyhc 

q,n,v 



'p( n -<l) e - pH l <p{n-q) 



Z ((V-q)/V,n-q,V-q). 



(3.14) 



Remark 3.3 For s, a G 



q) e -sH l >p(n-q) 



pin-q) e -saW <p{n-q) 



We expand 



trace 



q,n,V 



trace 



1,n,V 



in a Dyson series in powers of P. If m > 1, the m th term is 



X m := /3 m / dsi / <is 2 . . . / ds m trace ^hc 

'o Jo Jo q ' n ' v 



e -/3H(l- Sl ) pe -f3H( Sl -s 2 ) p _ . 



. Pg—/3H(s m ~i-s m ) p^-/3Hsmjj 



Recall that P := V* c Pi q) V* c where 



P (?) 



n. v ' n 
Pv®/®---®/ + /®Pv®/ 



(g>i"(g> Py 



(3.15) 



has g terms, so in the above trace we have m instances of this form. Let P 



(9) 



P 



r place 



/, and let P r = V^P^V^ 



Then we can write 
• i 



o Jo 



X m = (T jdsx jds 2 ... I ds m • ■ • E 

n=l r m =l 







x trace oyhc 

1,n,V 



-PH(l-S!)p -f3H(s!-s 2 )p ...p p -PH(sm-i-s m ) p -j3Hs m jT 



In terms of the basis of 7~Cg C n v we may write the expression for X m as 



"q,n,V 



x m = r d Sl Us 2 ... d Sm - E E 

ri=l r m =l k ,... ,k m i°~k° i m ~k m 

kV^ 1 "' ^ 1 ; k 1 )(i 1 ; kV^ {si ~ S2) ^r 2 |i 2 ; k 2 > • • ■ 
• • • (i" 1 " 1 ; k m - 1 | e -^ (Sm - 1 - Sm) P r Ji m ; k m )(i m ; k m |e-^ s ™?7 ? |i ; k°) (3.16) 

where it is understood that the i summations are over Ay and the k summations are over 
A?+ q) - Note that for i ~ k 



P r |i;k) = — Y ,ir,l,--- ,ig); k) (3.17) 

i=i...y^ 
i^k; l^i\...i r ...i q 

where again the hat symbol implies that the term is removed from the sequence. 



Consider one of the inner products in the expression ( 13. 16ft for X m , using ( 13.111) and (13.171) 
above. For i ~ k and j ~ k': 



(i ; k|e-^P,|j;k') = -— J2 (i\Ui,---S,l,---,j q ))(k\e-e sH '\k>) 



V 

1=1.. .v 
l$k';l^ji,...j r ,...j q 



S ilh ...S irjr 5 irl ...^(kle-^lk')- 



V 

1=1.. .v^ 
l<£U';l^jl,...Sr,-,jq 



In summing over I we replace / by i r and the result is non-zero only if i r ^ k' and i r ^ 
ji, . . . , j r , . . . ,j q . However this last condition is not necessary because if i r = j s (s ^ r) then 
j s 7^ i s and we get zero. Also if for some s ^ r, i s e k' then once again j s ^ i s . We can 
therefore replace the condition i r ^ k' by i ~ k'. Using I for the indicator function, we have 



e -Pqs 



(i;k|e-^P r |j;k') = — ^...i,, . . . 8 iqjq (k\e^ sH \k')T^ w) 



V 

Now if we sum over j ~ k', with i ~ k and for a fixed r: 



5- 5- 8- /klP (n " 9) e _/3sHi P (n " 9) lk'\ 



^(i;k|e-^ s P r |j;k')(j;k'| = ^—^vt^^^t^) 



X 



j r =l...V 

>£k'Ui\{>} 



tb 



It is convenient to define the operation [r, x](i) which inserts the value of x in the r 
position of i instead of i r . So for example taking the ordered triplet i = (5,4,1), then 
[2, 8] (i) = (5,8,1). For brevity we shall denote the composition of these operators as 
[r k ,x k ; . . . ;r 2 ,x 2 ;r 1 ,x 1 \ := [r k ,x k ] o • • • o [r 2 ,x 2 ] ° [ri,Xi]. 

Thus the final term in the above expression may be rewritten as ([r, j r ](i); k'|. 



Performing two summations for fixed r\ and r 2 we get: 



ii^k 1 i 2 ~k 2 

e -Pq(s+t) 



E E ^ k °i^ e "^ i0 ^°' kl ) 



^ 2 



-/3g(s+t) 



V^2 ^ 



(k°| P i0 e^" V i0 Ik 1 ) 



x (k 1 ! P [rii4i](i0) e -^ [ri4lI(i0) V [ri , ihm (^(h.^jr^g^jk 2 ! 

due to the fact that j^k 1 ) = if i\ x G k 1 . We may apply this to all inner product 

terms of (13.161) except the final one. Note we sum over the V sites of the lattice, with certain 
points excluded in each case. 

For the final inner product of (13. 16[) we obtain: 
(K,C; ••• S^^^r^zJJ^sk-le-^^l^k ) 

= e-^([r m ,C;...;^ 

= e-^(k-|e-^ [rm ^ ; ''- ;r2 ^ ;ri4lKi0, |k°) ([r m ,C; ... ; r 2 , z 2 ,; r 1; ^(i )^ ) 



e (k |/> m ,c m ;...;r- 2 ^ 2 ;r 1 ,ii i ](i0)e Hrm,i^i...;ra,^ 3 ;ri,il 1 ](P)|K ) 

x ([r m ,C; ••• ;r 3) i?,;r 1) tJ 1 ](i )|Z7 ff i >. 



Applying this to the whole tracial expression of (13.161) we obtain 



am 
ym 



X m = e-^ E E E 



E 



E 



k°...k™ i° i\0 a \{i a ri } i? 2 ^[ri,ii 1 ](i°)\{ii 2 } i™J[r m _ l! C 1 1 ;...;r2^;ri 1 «y(i»)\{C m 1 } 

(k^^ioe-^-^^Violk 1 ) 
(k 1 |P [ri4i](i o ) e-«---) 5[ri ' 4lI(i °V [riA](i o ) |k 2 ) 

\ k \> ;ri,ii ](iO)e /^ 2)i 2 ;ri)i i ](i 0)|k ) 



\ K K['"m 1 j^ r! ;-;r2 1 j? 2 ;r-i,ii 1 ](i ) e '[?•-' 
([r m ,C m ; ...;r 2) ^;r 1 ,^ 1 ](i°)|C/ 9 i ) 



,iP 2 ;ri,ii 1 ](i°)l k °) 



-/3g 



I" 



E 



i° i^U^} ^[n,iy(i°)\« 2 } l ™J[r m _ 1 ,C ra -_ 1 1 ;...ir 2 ,«? 2 ;ri, 1 i 1 ](i°)\{C 1 l 



trRCG fyhc 

^(71-4),^,+ 



(i°) e 



(i°) 



■Pi 



r-m,^;- ;r- 2 ,i2 2 ;ri,ii 1 ](i°) e 



P 



r-m,^;- ;r 2 ,i2 2 ;ri,ii 1 ](i°) 



From the Holder inequality (see Manjegani [ID]), for finite dimensional non-negative matrices 
Ai, A 2 , ■ ■ ■ , ^4 m +i we have the inequality 



m+l 



[trace (Ait4 2 . . . A m+ i) | < trace | Ax A 2 ■ ■ ■ A m+ i\ < J [ (trace A p k k ) 



j_ 

Pk 



k=l 



where YI=i ^ = 1, Pi > °- 



Set pi = p 2 = ^ 

trace 



S1-S2 



, • • • , Pm 



— j Pm+i = — ■ Taking the modulus of the above 



— 1 b m 



tr£lCG -xyhc 



p i oe-^ i ^ 1 --)p i oP [n ,, i](i o ) e-^ [ri ' i - ](1) (---)p [ni4i](i o ) 



< trace W hc 

rL (n-q),V,+ 



' ' trace oyhc 

rl (n-q),V,+ 



Pi 



~\ l-si 



P[r 1 ,iiJ(iO)e 



Pi 



n,4J(i°) 



)5[rm,iJ!L ; — ;r 2 ,j^ ;r!,i r ](i°) 



' [rm,i™ m i ... ;r2,«2 2 ;n,ii 1 ](i ) t; ' [r m ,i™ ; ... ;r 2 ,i^ 2 jn.i^ ](i°) 



try q * 



Since the trace is independent of the V — q sites [7*1, ](i°), 
; r 2 , z' 2 . ; r±, z* ](i )}, and therefore using Remark 13.31 the product of all the trace terms above 
is equal to 



trace oyhc 

n -(n-q),V,+ 



with 1 = {V — q+ 1,V — q + 2,. . . , V} and from Remark 13.21 



trace -^h 



(n-q),V,+ 



P.e-^P, 



Z ((V-q)/V,n-q,V-q). 



(3.18) 



Consider the sum 

EE E 

([^C;---;^V 2 2 ;ri,4j(i )|f/ g i°). (3.19) 

^[r ra _ 1 ,i^_ 1 1 j...;r a ,i5| a ^i I ii 1 ](i'>)\{Cr 1 } 



If {ri,r 2 , . . .,r m } ^ {1,2, . . . ,q}, then |[r m ,i"; . . . ; r 2 , i 2 r2 ; r x , ijj(i )) is of the form 

J2> • • • ijnii *ni+l' ' ' ' ' ^ri2'^™2+l5 • • • ) Jri3) ? 'ri3+l> " " " ' ^ri4'i".4+l) ) 

where {rii,n 2 , ■ ■ ■ } is a non-empty ordered set of distinct integers between and q. This 
state is clearly orthogonal to U q i° for any q. Note that this situation does not arise if q = 1. 
Note also that this is always the case if m < q. 

We may bound the remaining sum corresponding to terms for which {r 1; r 2 , . . . , r m } = 
{1,2,. ..,<?} by 

^E E E E ([r m ,C;...;r 2) .? a ;r 1) <y(iW>- 

i° ^=1 ia,=l ij^=l 
v v ' 

where [r m , ; - ■ ■ ; T"2, ^ 2 ; J"i , « r ] 
has distinct indices 

Observe that in this case | [r m , z^j ... ; r 2 , z'J: 2 ; ri, z*J (i )) is independent of i° so we may take 
it to be 

|[r m ,C; ... ;r 2 ,il 2 ;r 1 ,i 1 ri ](s )) 

where s° = (1, 2, 3, . . . , q). Then we can interchange the i° summation with the others, and 
for each choice of z* , if 2 , . . . , i™ m there exists only one possible i° G such that 

([r m ,C; ... ;r 2) ^;r 1) <J 1 ](i )|l7 fl i >^0 

So we may conclude that 



EE E 

i° i 2 r 2 <t[ri,i 1 n )(i )\{i 1 r 2 } 

<[r-m,C;---;^V 2 2 ^i>4](i )l^°}<^ m (3-20) 



Applying this, we see that the modulus of the integrated m th term of the Dyson series may 
bounded above by 



am i 



ml V 



ri=l r m =l 



EE E 



E 



([w™; ■■■ ;^4;n,i'](i )|f/,,i ) 



Vm Fl r m-lA 



m — 1 ' > r 2i V 2 >' il^riJV 1 ; \ \"r m 



am i 



ml V 



ri=l r m =l 



e ((y_ g )/v;n- g ,y- g ). 
m! 

Noting that the zeroth term of the Dyson series is 

Xn = trace W hc 

n q,n,V 

we may re-sum the series to obtain 



' tr£lCG n-/hc 



<e-toZ p {(V-q)/V,n-q,V-q) £ 



m=l 



ml 



Thus 



te(q)-<$(q)\ 



V 



1,n,V 



U q e-^,v 



tr£lCG ovhc 



U q e~P H 



Z (n,V) 



I" 

— 1 



m! 



e^ - 1) 



Zf,(n,V) 

' v ' m=l 

^((V-g)/V,n-g,y-g) 



3.3 Proof of Proposition EO 



Recall that we have 

Zft{n — q,V — q) = trace H hc v + [ e ~ /3Hn ~i' v -i} 



while 



-0(n-q) 



trace W hc 



n — q,V — g,+ 



n—q V—q rt — q 



(3.21 



^ n-q,V-q) = trace W hc_ y _ [e ^(v-d/v.n-ftv-a] 



-/9(n-9) 



trace 



n — q,V — q,+ 



o/ V — q \-phc p^ — Q'phc 



(3.22) 



(3.23) 



Comparison of (13. 21ft and (13.221) yields 

Z, ({V-q)/V, n-q, V-q) = e^v>-<*) Z^ijx - q,V - q) 

and thus we have to analyse the following ratio: 

e-^M^ (ts) (n - q, V - q) 
Z (n,V) ' 
Penrose in [T3] gave an explicit expression for Z (n, V): 

min(n,V—n) 

Z p (n,V)= Yl *(T,n,V,P), 

r=0 

where 

z{r,n,V,p) := Oy^^Yj LJ exp {~y ~ r ' + r + n ' ~ ™] } • 

He also proved that if h v : [0, min(p, 1 — p)] — > R converges uniformly in [0, min(p, 1 — p)] as 
V — > oo to a continuous function h : [0, min(p, 1 — p)] — > R, then 



min(n,V-n) (Hp)i if P ^ Pfi], 

hm — -— ^ Mf)«(r,n,V,/3) = ^%^), if p € [fp, 1 - (3.24) 



We wish to express the ratio in (13.231) in the form of the lefthand side of ( 13.241) . We have 

min(n-q,V-n) /y_ _ 2 r + \\ fV — \ 

Z^(n-q,V-q)= g ( v-q-r + l ){ r J 

x exp / — ^ [r(V — q) — r 2 + r + (n — q) 2 — (n — q)~\ \ 

For the case p > |, for large V, n — q > V — n we must sum from zero to V — n and a 
straightforward calculation then gives 



V-n 



v 

r=0 



e-M^Z^^n -q,V-q) = Yl M£) z ( r > n > V > P) 

where 

,'l-2x- (q-l)/V\ ( 1-x + l/V 

h v (x) 



l-2x+l/V J \1 - x - (q - 1)/V 

q-1 

X 

s=0 



n ( 1 1 - s /v v ) exp {(3q [x+p ~ i/v]} ■ (3 - 25) 



Therefore 



h(x) = lim hy{x) = (1 - x) 9 e^ 2 ^. (3.26) 



It is clear that the convergence is uniform since h v (x) is a product of terms each of which 
converges uniformly on [0, 1 — p] for p > |. Thus 

e-m^Z^^n-q.V -q) ( (1 _ Pft )< e tfQ>?+ P ) if p £ (1/2, 1 — pg\, 
hm 7 - = < 

^(n, 1/) [pete if p e (1 - p„ 1]. 

Note that using the relation 



l-2pp V P/3 
we get 

(1 - p p ) q e qf3(p e +p) = pi e ^ 1+ P-P^ 

and therefore we have proved Proposition 13.31 for p > |. For the case p < | we have that 
ri — g < 1/ — n, the sum for q-Pv^-i) Z ^v-z^n — q,V — g) is up to n — q, and therefore we 
need to shift the index by q to get it into the required form. After shifting we get 

.-**<"-% V) (» - «, k - ,) = £>, », v, « ( Vl^+l 1 ) 

r=q ^ ^ 

x K^-l)(^-2)---(r-g + l) (f3q lv } 

X nV-l)(V-2)---(V-g + l) eXP \V rJ 

Note that summand is zero if we put r = 0, . . . ,q — 1. Thus we may sum from zero to n to 
get as before 

n 

e-rf^Zp^ (n-q,V-q) = J2 Mf) *fa "> V > P) 

where this time 



V 

r=0 



, , , fl-2x+(q + l)/V\ 14 fx-s/V\ rn . „ , , 



so that 

h(x) — lim = x 9 exp{/3g(l + p — x)}. 



n,V — »oo 

n/v=p 



Convergence is again uniform on [0, p] for p < \ and therefore 

e-^-% (V) (n- g ,T/- g ) f^eftr ]£pe[O,p ), 

hm t r = < 

^(n, 10 \p« e ^( 1+ "-^) if p e \p p , 1/2), 

proving Proposition 13.31 for p < ~. The case p — ~ is more delicate because the first term in 
(I3.27P does not converge uniformly. We can write (taking V = 2n) 

~ q ~ 

h 2n (r/2n) = h 2n (r/2n) + — ——h 2n (r/2n) 

2[n — r) + 1 

where 

9-1 



h2n(x) = J! (l3T7§0 ex P V* P/2 - *]} ■ 



Clearly h2 n {x) converges uniformly on [0, 1/2] and therefore 

1 n ~ 

We thus have to show that 

1 h,J^ 

hm 



n^oo Z /3 (n, 2n) ^-^ 2(n — r) + 1 
Since h 2n (x) is bounded, by C say, 

lim — - V z(r,n,2n,(3) < lim = 0. 

n-oo Z g (n, 2n) ^ 2{n - r) + 1 v ' ' ' M; ~ n-°o 2nV4 

r<n— n 1 / 4 

On the other hand one can prove that for n — 2n 1 ^ 2 < r < n — rt}l 2 and r' > n — n 1//4 

In z(r, n, 2n, (3) — In z(r', n, 2n, (3) > - In n 

8 

for n large, so that z(r' ,n,2n, f3)/z(r,n,2n, (3) < 1. Therefore 

^2z(r,n,2n,f3) 

1 h 2n (j^) r > n - n i/4. 

> — =f= z[r,n, 2n, p) < lim C 



lim 

«-ooZ/n,2n)^ i/4 2(n-r) + l n^oo ^z{r, n, 2n, /?) 

n—2n 1 / 2 <r<n—n 1 / 2 

maxz(r, n, 2n, (3) 

< lim ——- — r ~ n — - < lim —r- = 0. 

mm n 1 ' 4 mm z(r, n, 2n, (3) ^m n 1 ' 4 

n— 2ra 1 / 2 <r<ra— n 1 / 2 



3.4 Proof of Proposition 13.41 



Recall that 



Considering the trace over TL q c n v , expanding it in terms of its basis {|i; k)} and using Remark 
13. II above, where i ~ k 



trace U^KM = ^ ^<i; k\U q e~^ n 



hc^_p(™-«)jphc 



i;k) 



k i~k 



k i~k 



k i~k 



For q > 1, an element of the basis of the unsymmetrised g-space may be written as an 
ordered g-tuple i = (i 1 ,i 2 , . . . , i q ) where the ij's are all distinct. Then we may write 



(U q i\i) = (UJei^ei 



e k ® e it | e it <g> e i2 <§ 



e «,) 

= o. 



Hence c"(g) is non-zero only if q = 1. 

For the second statement, note that we may re-express c"(l) as follows: 



1 1 



;trace 



-P i 



v 



Z 9 (n,V)V^ 



^^(k|e^|k) 



-p 1 



E trace W hc 
rl n-l,V,+ 



Z,(n,V)V^ 

_ Z,(^,n-l,V-l, 
" 6 Z p {n,V) 

and the result follows from Proposition 13.31 
4 ODLRO 



The one-body reduced density matrix for x, x' £ A v may be defined as 

1 



D Pt n,v(x, x ) :— (a x a x ,) - 

where for £ H v , K x>x >4> = (e x >\(j))e x . 
Penrose showed that for i/f, 



— — r-r-trace oyhc 

Z g (n,V) n "> v -+ 



lim D^ njV (x,x') = p% 

n,V — »oo 

n/v=p 



(4-1) 



that is, whenever Bose-Einstein condensation occurs, there is Off-diagonal long-range order 
as defined by Yang [20J. It has been argued and proved in some cases (see for example [18] 
and [5]) that in the expansion of Dp n v (x,x') in terms of permutation cycles, only infinite 
cycles contribute to long-range order. Here we are able to show this explicitly. 

By the proposition in Appendix [A], we have 

n 

Dp >ri7V (x,x') = K XtX >) 

9=1 



where 



C?(q;K x 



1 



— — — r trace -ryhc 

Zf,(n,V) n ^ v 



(4.2) 



{K XiX , <g> I <g> I <g> . . . <g> I)U g e t* H «v 

Note that this is equivalent to the expansion of o~ p {x) in equations (2.14) and (2.16) in [18] . 
Applying the argument in Subsections 13.21 and 13. 3^ we can show that 



lim C^{q-K x y)= lim C%(q;K x 

n,V — ► oo n,V —^oc 

n/v=p n/v=p 



(4.3) 



where we take 

C^{q; K X)X ,) = / trace W hc [(^^ ® I ® J ® . . . ® I)U q er pH ^ 
The only difference is that instead of equation (13.191) . we obtain 

EE E - E 

([r m ,i£; ... ;r 2 ,v 2 2 ;r 1 ,^ 1 ](i°)|(^®/®---(8)J)C/ g i ) (4.4) 

whose treatment is similar but slightly more complicated, as detailed below. 

Let q > 1 and consider the case {r±, r 2 , . . . , r m } 7^ {1, 2, . . . , q}. When 1 ^ {r±,r 2 , . . . , r m } 
we obtain inner products of the form: 

(i\\K x ^il){j 2 ,j^ ■ ■ ■ , .7 ff |i3,i°, • • • 

where jk 7^ i\ for all k by the hard-core condition, implying the second term is zero as 
j q 7^ i\. On the other hand, when 1 G {ri, r 2 , . . . , r m }, then there exists at least one 
/ ^ {7*1, r 2 , . . . , r m }, yielding an inner product of the form 

(ji\K XjX d° 2 )(j 2 , . . . . -Jqli^il, ■ ■ 

which also results in the second term being zero as = 0. Note that the above cases 

do not occur for q = 1. 

For the case {7*1, r 2 , . . . , r m } = {1, . . . , q}, as before, the remaining sum may be bounded by 
a similar expression whose summations have slightly relaxed restrictions. Also the left hand 
side of the inner product is independent of i°, so again denoting s° = (1, 2, 3, . . . , q), we have 

V V V 

ra< E E • • E E<[ r -C;---;^ 2 ;ri,^ 



L r 1 — 1 ?" 2 — 1 1 



where [r m , ; ■ ■ ■ ; T i > ; r i , ] ) 

has distinct indices 



and as there is only one possible value for each i®, i®, i®, . . . , i q giving a non-zero summand, 
we can bound above by 

< E E ■•• E E«J^.^5) = v m - 1 E XXi*--© = vm - 1 

i\=l i= 2 =l i™ m =l;o=i i* fc =l ig=l 

where G [1, m] is the smallest number such that = 1, and for any x,x' G A v . Thus the 
entire sum (14.41) is bounded above by V m ~ l . Therefore one can conclude the argument of 
Subsection 13. 3\ proving (I4.3p . 

Moreover, following the reasoning in Subsection 13.41 we can then check that for q > 1 and 
x 7^ x', Cy(q; K XjX r) = 0, since for q — 1, (ei\K XtX rei) = 0, and for q > 1 

{(K X}X , ® I ® • • • ® = ((i^e^) <g> e i3 ® • • • <g> e i? ® e ix | ® e i2 ® • • • ® e iq ) = 



as the ii's are all distinct. So we have that 



lim C%(q]K XtX ,) = 

71,, V — »00 

n/v=p 



and that 



lim lim V" C^{q;K x ^ x ,)= lim D 0>n>v (x, x') = p c . 

O — >00 n,V — >oo n,V — ►oo 

n/v=p 



>00 n,V — >oo 

n/v=p q=Q+i 



A Appendix: Expectation of Operator in Terms of Cycle Lengths 

Proposition A.l Given an operator A on 7i v , the expectation of A may be expressed in 
terms of cycle lengths 

n 

{A (n) ) = J2c-(q;A) (A.l) 

9=0 

where 



^ A )=Z^V) traCe Kn.V 



(A® I® I)U q e 



-mi 



(A.2) 



Note that c"(g) = Cy(q; I)/V where c"(g) is as defined in Proposition 13.11 
Proof: 



(A {n) ) = — —-—trace H hc 



Z p (n, V) n\ 



^trace^ \a^U^ h ™ 



using the facts that [U n ,H^ c v ] = 0, [f/ 7r ,'P^ c ] = and by the cyclicity of the trace. Note we 
can simplify this expression by the following method: 



trace avhc 

"n.V 



trace 



he 

n,V 



*=1 ith position 

n 

^ f/(ii)(A ® J ® • • • ® J)t/ (li) ^e-^V 



i=l 



where Un^ represents the transposition (1 i), so using cyclicity of the trace again 



trace 
trace H 



he 

n.V 



^2(A ® I ® • • • ® I)U {li )U n e P H % v U m 



i=i 



he 

n.V 



® J ® • • • ® I)U {li) U^U {li) e 



= n trace W hc^ | (A ® J 
where 7r' = (1 i) ir (i 1) using (13.101) . Thus 



(A (n) ) 



(n-1)! 



trace W h C (^4 ® J ® • • • ® J) C^e 



(A.3) 



Given distinct indices 12, ■ ■ ■ ,i q , let 

£'(12,13, • • -h) = I 71 " e s n ■ 7r(i m ) = i m+ i, 1 < m < q with z a = w(i q ) = 1 j. 
Then for any 7r £ S 1 '^, 23, • ■ ■ iq), there exists a 7r' £ S n - q so that one can write 



trace 



(A ® I ® • • • ® I)U n e 



trace 



he 

n.V 



{A® I ®---®I){U q ®U^)e 



-/3H 



The set S^(i 2 , 13, ■ ■ ■ i q ) form a partition of the set of permutations where 1 belongs to a cycle 
of length q. There are such sets. Then 



{n-q)\ 



(AW) 



Z (n,V) (n- 
1 1 



1)! E 



trace ovhc 

n n,V 



tt€S„ 
n 



(A ® I ® • • • ® I)U^ pH ^v 

(A ® J ® • • • ® J)(Z7 ff ® U*)e- pH %v 



Z,{n,V) {n-l)\^ {n-q)\ 

1 A 1 ^ 



1 V 



trace ^ c(H (g)^(^« ?)) 



{A® I 



I)U q e 



-m. 



n,V 



and recall that H h „ c n v := Vl\nf ® W 



#(n-g)> 



□ 
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